Pattern of Reaction Diffusion Front in Laminar Flows 
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Autocatalytic reaction between reacted and unreacted species may propagate as solitary waves, 
namely at a constant front velocity and with a stationary concentration profile, resulting from a 
balance between molecular diffusion and chemical reaction. The effect of advective flow on the 
autocatalytic reaction between iodate and arsenous acid in cylindrical tubes and Hele-Shaw cells 
is analyzed experimentally and numerically using lattice BGK simulations. We do observe the 
existence of solitary waves with concentration profiles exhibiting a cusp and we delineate the eikonal 
and mixing regimes recently predicted. 



The motion of interfaces and the propagation of fronts 
resulting from chemical reactions occur in a number of 
different areas Q , including population dynamics 0, 
and flame propagation Q. It is known that autocat- 
alytic reaction fronts between two reacting species prop- 
agate as solitary waves, namely at a constant front ve- 
locity and with a stationary concentration profile 0, . 
The important issue of the selection of the front velocity 
was addressed earlier on, but only a few cases are well 
understood, such as the pioneering works of Fisher 
and Kolmogorov-Petrovskii-Piskunov 1^ on a reaction- 
diffusion equation with second-order kinetics 0, 0, 0- 
The effect of advective flow (inviscid and/or turbulent) 
on reacting systems was analyzed extensively in the prop- 
agation of flames in the context of combustion |j, 
On the other hand, advective effects on the behavior of 
autocatalytic fronts have been only recently addressed 
0,110, [nf. B. F. Edwards studied theoretically the 
effect of a 2D laminar flow on an autocatalytic reaction 
front between two infinite planes separated by a gap b. 
In this geometry, the velocity profile is unidirectional in 
the direction z of the flow and is given by Poiseuille's 
equation, if = /7m(1 — C^)^ where Um = 1-5 U is the 
maximum velocity, U is the mean velocity, ( = 2x/bis the 
transverse normalized coordinate and ^ is the unit vec- 
tor parallel to the flow, chosen as the direction of the front 
propagation in the absence of flow (see below). Consider 
the iodate-arsenous acid reaction described by a third- 
order autocatalytic reaction kinetics 0, 0, : 



dc_ 
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Tf -^0 = DmAC + aC\l - C) 



(1) 



where C is the concentration of the (autocatalytic) re- 
actant iodide, normalized by the initial concentration of 
iodate. Dm is the molecular diffusion coefficient, and a 
is the reaction rate kinetic coefficient. In the absence 
of hydrodynamics (if = "O^), Eq^ admits a well-known 
solitary wave solution with front velocity Vq = ^/aDm/2 
and front width Lq = Dm/Vo 5, 6]. The use of these two 
quantities to normalize velocities and lengths in EqQ 
leads to two independent parameters rj = b/2Lo and 
s = U /Vq. Reference llj investigated numerically the 
solitary wave solution of EqQ and particularly its nor- 
malized front velocity, v = Vf/Vo^ as a function of £, 



for different values of 77. Of interest are the following 
asymptotic predictions: 

In the narrow-gap regime (77 ^ or £ ^ 0), it was 
found that v = 1 -\- e. Namely, when Lq » 6, mixing 
across the gap is significant, the concentration front is fiat 
and advected by the mean fiow, yielding: Vf = Vo -\-U. 

On the other hand, in the wide-gap regime (77 ^ 1), 
the front is thin and curved across the gap, and Eq^can 
be replaced by the eikonal equation: 



Dr, 



(2) 



where it is the unit vector normal to the thin front (ori- 
ented from reacted to unreacted species) and tz the front 
curvature. In this regime, to leading order and neglect- 
ing the local front curvature, the front velocity is given 
by the simplified ID eikonal equation: 



Vf = Vo/cosO^U{C) 



(3) 



where is the angle between it and the fiow direction and 
/7(C) is the advection velocity. Under these conditions, 
reference ^ predicted two behaviors depending on the 
fiow direction: For a supporting fiow (e > 0), Vf = + 
Um^ which means that the front is advected at the largest 
possible velocity. The front shape across the gap is then 
given by the solution of Eq|31 For an adverse fiow (e < 0), 
Vf = Vb, which also represents the maximum algebraic 
velocity one could have expected physically. The front 
is perpendicular to the walls (6> = at = il), and 
presents a cusp in the middle of the gap (discontinuity of 
6^ at C = 0). Here, the adverse fiow elongates the front but 
does not slow it down. Note that curvature effects {Dmi^ 
in Eq|2|) smooth the cusp, but do not otherwise alter 
these predictions. Note also that similar features would 
occur for other kinetics such as FKPP or Arrhenius ones 

The objective of the present letter is to experimentally 
test the above 2D predictions using two different devices, 
namely Hele-Shaw cells and cylindrical tubes. The case 
of the Hele-Shaw cell, consisting of two parallel plates 
separated by a gap b small compared to the other di- 
mensions, is supposed to be quantitatively addressed by 
Alternatively, the case of the cylindrical tube of 
inner radius a (in which the fiow field is also described 
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by Poiseuille's equation, with Um = 2[/ and ( = r/a) 
represents a genuine (axisymmetric) 2D situation. Ex- 
periments in Hele-Shaw cells are discussed with the help 
of lattice BGK simulations [i3 of EqlUfor a 31) flow. 
In the experiment, the front is detected by using starch, 
at small concentrations, which reacts in the presence 
of iodine leading to a dark blue signature of the front 
passage 

11113. 

First, we consider the reaction in the 
absence of advection by the flow (if = "O^). As ex- 
pected, we do observe solitary fronts propagating with 
flat shapes. In accordance with 5, 6], their velocity is 
1^ ~ 0.02 mm/ from which one can estimate their front 
thickness Lq = Dm/VQ ^ 0.1 mm {Dm ^ 2.10"^ m^ /s). 
Because the reaction products have a lower density than 
the unreacted species, the hydrodynamically stable situ- 
ation corresponds to descending fronts in vertical tubes. 
In the following, we focus on the interplay between ad- 
vection and propagating fronts. To minimize the effect of 
density contrast, we studied the propagation of buoyantly 
stable fronts in small cells. We used vertical Hele-Shaw 
cells of size 6 x W = 0.1 x 1, 0.2 x 4, 0.4 x 8, 1 x 15 mm^ 
and circular capillary tubes of radius a = 0.3, 0.58, 0.88 
and 1.9 mm. A constant advecting flow, upwards or 
downwards, was fixed by a syringe. Note that these cells 
are small enough to prevent flattening of the front due 
to buoyancy, but large enough to enable a constant flow 
rate injection with our injection device. The average ve- 
locity of the imposed flow ranged between and ~ 60 Vq. 
In a 3D Hele-Shaw cell, the flow velocity profile is uni- 
directional and depends on the two transverse coordi- 
nates, X and y [l^- The profile across the gap is al- 
most parabolic with a gap average value uniform over the 
width except in a boundary layer of order 6, within 
which the velocity vanishes (see the gap-average profile 
on Fig^. For the three aspect ratios studied, W/b = 
10, 15, 20, we can estimate from ^,Um/U = 1.60, 1.57 
and 1.55, respectively. We observed solitary waves in the 
whole range of flow rates investigated. Typical fronts are 
shown in Fig^in the plane of Hele-Shaw cells (top dia- 
gram) and in tubes (bottom diagram). For each geom- 
etry, two adverse flows (on the left) and two supportive 
flows (on the right) are displayed. The front shape al- 
ways points toward the same direction as the underlying 
flow field, while its distortion increases with flow inten- 
sity (recall that the fronts are flat in the absence of flow, 
U = 0). The triangular shapes observed, in the case 
of adverse flows, are reminiscent of premixed flames |T5| . 
The two graphs in Fig^show the normalized front veloc- 
ities V = Vf/Vq versus e = U /Vq, measured for different 
sizes of the Hele-Shaw cells and tubes. For each data set 
(given 77), the front velocity increases linearly with the 
flow rate, but with a slope different for supportive and 
adverse flows. This change of slope is in accordance with 
im as well as the observed evolution of the slopes with 
77. In addition, most of the data fall in the domain de- 
limited by the asymptotic regimes described above. The 
exceptions correspond to supportive flows in the smaller 
tubes (1 < 77 < 6), which, unlike Edward's 2D numerical 
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Figure 1: Normalized front velocity v versus normalized flow 
velocity s (s < 0: adverse flow, s > 0: supportive flow) for 
different normalized sizes 77. Top: Hele-Shaw cells of differ- 
ent normalized thicknesses and aspect ratios (77 = 6/2Lo, W/b): 
0(0.5, 10); #(1,20); 0(2, 20); ■(5,15). Bottom: Circular tubes 
of different normalized radii (ry = a/Lo) : ■(3); 0(5.8); •(8.8). 
The full and the dashed lines correspond respectively to the mix- 
ing regime (77 0) and to the eikonal regime (77 oc). Exper- 
imental pictures: From left to right s = -4.8, -2.4, +2.4, +4.8 
for Hele-Shaw cells and e = —6.7, —1.9, +1.9, +6.7 for tubes. 



data, fall below the mixing regime (predicted for 7^ <C 1). 
This difference needs to be further analyzed, given that 
it is difficult to achieve experimentally a very low con- 
stant flow rate injection in the small tubes. At the same 
time, some authors lOj have suggested that the mixing 
straight line should be higher. Using the Peclet number 
(Pe = Ua/Dm = £^), which compares the relative im- 
portance of advection and diffusion, and the Damkohler 
number {Da = aa/U = 2r]/e)^ which is the ratio of ad- 
vective to reactive time-scales, they predicted that 
for Da <C 1 (which is not attainable in our experiments), 
the front velocity should be the product of Vq by the Tay- 
lor dispersion factor 16], which accounts for the couphng 
between advection and transverse diffusive mixing. This 
factor would then enhance the front velocities. 
The front velocities measured in the Hele-Shaw cells are 
very close to the mixing regime. However, in the case of 
adverse flows, the measured values exhibit some depar- 
ture toward the eikonal regime {Vp = Vq) when either 
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Figure 2: Normalized front velocity v versus normalized flow 
velocity s for the smaller Hele-Shaw cell of size 1 x 0.1 mm^ 
(?7 = 0.5, W/b = 10). The full and the dashed lines correspond 
respectively to the mixing (77 0) and to the eikonal (77 00) 
regimes. 



77 or £ is increased. This trend is even more pronounced 
for the tubes, which present larger 77 values than for the 
Hele-Shaw cells, in accordance with predictions On 
the other hand, all the values measured for supportive 
flows, even large, {e up to 50 for 77 = 0.5 displayed in 
Fig|2|) fall on the asymptotic mixing regime predicted 
by the strictly 2D gap analysis {Vf = Vo -\- U for 
77 0). This is all the more surprising since the flatness 
of the fronts, expected in the mixing regime, was not ob- 
served in the Hele-Shaw cell plane (top right of Fig^. 
However, in our experiments, although the normalized 
gap T] = b/2Lo introduced in the 2D gap analysis |0 is 
small (0.5,1,2,5), the normalized width W/2Lo is large 
(5,20,40,75). Extrapolating Edwards's 2D gap analy- 
sis to our 3D case would suggest that our experiments 
combine a mixing regime across the gap with an eikonal 
regime across the width. Under these assumptions, the 
shape and velocity of the front would obey an equation 
similar to Eq|21 namely 

Vf = {Vo + Dmf^)/ cos 0^ U^^'iy) (4) 

where the effective advection velocity U'^^ {y) is the gap- 
averaged velocity defined in , and where k, represents 
now the curvature of the 2D front curve observed in the 
plane of the Hele-Shaw cell. The front velocity for e > 
would then be set by the maximum Ulf of the profile 
U'^^{y)^ found in the middle of the plane (for 6{y = 0) = 
and Hi{y = 0) = 0). The so-obtained maximum front 
velocity, Vf = Vq -\- UfP , expected in the asymptotic 
width-eikonal regime {W/2Lo » 1) can be compared 
to the velocity, Vf = + /7, expected in the width- 
mixing regime {W/2Lo << 1). One finds that these two 
asymptotic velocities would be equal in the Hele-Shaw 
limit {W/b — > 00), and are actually very similar in our 



experiments (as UjP /U = 1.07, 1.05 and 1.03 for the 
three aspect ratios used). This could justify that the 
parameter r] = b/2Lo introduced in the 2D gap analysis 
[llj . actually controls the front velocity in 3D Hele-Shaw 
cells. 




Figure 3: Calculated concentration fronts in the plane 
of the cell, obtained by 3D lattice BGK simulations 
(full lines) and by integration of the 2D eikonal Eq^] 
(dashed lines), for 77 = 1 and for one adverse flow (top: 
s = —4.8) and one supportive flow (bottom: s = 
4.8). The front velocities so-obtained are Vf/Vo = 
—3.3 and 6.16, respectively. The dotted-dashed line 
is the front obtained by the 2D eikonal, when the 
flow velocity is slightly modified (2%) to mimic the 
effect of the meniscus which appears on the top of the 
reacted mixture. 



We tested the ability of the full description and the 
simplified one, given respectively by Eqs. Q and ^ to 
account for both shape and velocity of the experimental 
fronts. As lattice BGK simulations have been used to ob- 
tain the solutions of Eq^ we have first validated this nu- 
merical method, by reproducing Edwards's results on 
the shape and velocity of the fronts propagating between 
two infinite planes {2D simulations). Then, 3D lattice 
BGK simulations of Eq^and numerical integration of the 
2D eikonal Eq^were performed, using respectively the 
analytical stationary 3D flow field given by 14], and its 
gap-average U'^^{y). The front obtained with the latter 
method was compared to the iso-concentration C = 0.5 
of the gap-averaged concentration map produced by the 
lattice BGK simulations. Fig |31 displays these fronts, for 
the same parameters (e, 7^, W/Lq) as in one typical exper- 
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iment, in the cases of adverse {e = —4.8) and supportive 
{e = +4.8) flows. Note that in these cases of interest, for 
which W/Lq is finite, the integration of Eq^Jrequires the 
value of the front velocity (thus fixing the value of the 
curvature at the integration starting point). Hence, 
EqEI is not fully predictive, but links the shape of the 
front to its velocity. The shape and velocity predicted 
by lattice BGK simulations and the ones given by Eq^ 
are found to compare fairly well with the experimental 
observations in the case of adverse fiow. However, for 
the supportive fiow case, although the front velocity is 
correctly predicted, the two numerical predictions, sim- 
ilar for the shape, fail to account for the experimental 
observations. We believe that this discrepancy might be 
due to an alteration of the fiow velocity profile caused by 
a triangular meniscus which appears on the top of the 
solution in our supportive fiow experiments. From its 
shape and its distance to the front (typically several tens 
of W) , one can infer that the meniscus could introduce a 
few percents of excess fiuid velocity in the middle of the 
cell plane. The resulting non-uniformity in U'^^{y) may 
account for the rounded shape observed in the supportive 
fiow experiments (see Fig^. 




We have also analyzed the dynamics of the shape for- 
mation in the case of adverse fiows. Fig^ displays the 
time development of the iso-concentration C = 0.5, ini- 
tially fiat, toward the stationary triangular shape. The 
sequence shows that an early determination of both the 
final front velocity and the final angle is achieved as 
soon as the profile is altered over a typical distance b from 
the side walls. This supports the contention that W plays 
no role in the determination of both shape {0) and ve- 
locity in the regimes under consideration. This was con- 
firmed by simulations in wider lattices which produced 
the same values of the velocity and 0. Thus rj = b/2Lo is 
effectively the relevant parameter in Hele-Shaw cells (for 
which W/b > 1). 

In conclusion, we have performed experiments and lattice 
BGK simulations of autocatalytic reaction fronts in lam- 
inar advective fiow fields in Hele-Shaw cells and circular 
tubes. Solitary waves were observed in the entire range 
of fiow rates. For fiows adverse to the chemical front 
propagation, we observed cusp-like fronts in tubes and 
triangular fronts in the plane of Hele-Shaw cells. Our 
measurements of the front velocity agree with the 2D 
asymptotic predictions in the limiting cases where 
either diffusion overcomes reaction (77 <C 1) or it is negli- 
gible (7^ ^ 1). It would be interesting to extend the range 
of the cell sizes. Larger cells could be used to study the 
buoyancy stabilizing effect and smaller cells within the 
scope of microfiuidics. 



Figure 4: Time evolution of the front from a flat shape to a 
triangular one in the plane of a Hele-Shaw cell. The fronts ob- 
tained with 3D lattice BGK simulations (lines) are compared to 
the experimental stationary front (left of the flgure). The flgure 
on the right displays the gap-averaged flow velocity proflle of 
the simulation. 
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